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Segregation by thermal diffusion of an intruder immersed in a sheared granular gas is analyzed 
from the (inelastic) Boltzmann equation. Segregation is induced by the presence of a temperature 
gradient orthogonal to the shear flow plane and parallel to gravity. We show that, like in analogous 
systems without shear, the segregation criterion yields a transition between upwards segregation and 
downwards segregation. The form of the phase diagrams is illustrated in detail showing that they 
depend sensitively on the value of gravity relative to the thermal gradient. Two specific situations are 
CN| , considered: i) absence of gravity, and ii) homogeneous temperature. We find that both mechanisms 

(upwards and downwards segregation) are stronger and more clearly separated when compared with 
segregation criteria in systems without shear. 
I— 2 t Keywords: transport processes/heat transfer (theory), binary mixtures, granular matter 
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I. INTRODUCTION 

• The understanding of the physical mechanisms involved in segregation phenomena is one of the most important 
challenges in the field of granular matter. Apart from its academic interest, the problem is of central interest mainly 
due to its practical relevance in many industrial processes (powder metallurgy, pharmaceutical pills, glass and paint 
I ' industries, • • • ). In some cases, it is a desired and useful effect to separate grains of different types (e.g., the separation 
of mined ores) , while in other situations the resulting non-uniformity is an undesirable property that can be difficult 
to control. However, in spite of its practical importance, the problem is not completely understood yet. This fact has 
motivated the development of fundamental theories that provide accurate segregation criteria in the bulk region of 
. the sample [l|. 

One of the most familiar phenomena concerning segregation is the so-called Brazil-nut effect (BNE): when a binary 
I , mixture composed by one large ball and a number of smaller ones is vertically agitated, usually the intruder (large 
particle) tends to climb to the top of the sample against gravity On the other hand, a series of experimental 

works [y, 0] have also observed the reverse buoyancy effect, namely, under certain conditions the intruder can also sink 
to the bottom of the granular bed. This effect is known as the reverse Brazil-nut effect (RBNE). Although several 
mechanisms have been proposed to explain the BNE/RBNE transition [3-fToj] . the problem is still open. Among the 
different competing mechanisms, thermal diffusion becomes the most relevant one when the sample of grains resembles 
a granular gas (for example, at large shaking amplitude). In this regime, binary collisions prevail and kinetic theory 
can be a quite useful tool to analyze granular systems. Thermal diffusion (or thcrmophoresis) in dilute [lll - [l3j and 
dense |14H19j granular mixtures has been a subject of current interest in the past few years. 

A granular gas in rapid flow regime can be achieved by shearing from the boundaries [20| . Shearing is in fact as 
£SJ ■ common as shaking in experiments with granular systems. Thus, segregation criteria for sheared systems are also 
of interest. Furthermore, experimental works in annular Couette cells j2lM24| have shown that granular materials 
■ segregate by particle size when subjected to shear. Nevertheless, in spite of the relevance of the problem, much less is 
known on the theoretical description of segregation in sheared granular systems. In effect, to the best of our knowledge, 
t-H . previous theoretical studies [25[ on the subject for dense systems have been based on a Chapman-Enskog expansion 
around Maxwellian distributions at the same temperature for each species (26j . But the use of these distributions 
. , can only be considered as acceptable for nearly elastic particles where the assumption of the equipartition of energy 
still holds. Moreover, according to this level of approximation, effects of inelastic collisions appear only through a 
J-j \ sink term in the energy balance equation and for this reason the expressions of the Navicr-Stokcs (NS) transport 
- - ■ coefficients for the mixture are the same as those obtained for ordinary gases (elastic collisions). On the other hand, 
the use of the NS description to analyze segregation in steady granular flows is a serious limitation since the NS theory 
heavily fails beyond the quasielastic limit [271 ]. In addition, this is specially true in flows where viscous heating is 
exactly balanced by inelastic cooling, for which the granular flow is inherently non- Newtonian [28|. For this kind of 
flows, there is a special case of null temperature gradient, called simple or uniform shear flow (USF). This flow has 
received a great deal of attention in the past years and is the reference case study for granular flows (20j . 
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The aim of this paper is to analyze segregation by thermal diffusion in a binary granular mixture under USF in 
the framework of the inelastic Boltzmann equation. Due to the complexity of the general problem, here we consider 
the special case in which one of the components is present in tracer concentration. The tracer problem is more 
amenable to analytical treatment since there are fewer parameters than in a binary system. At a kinetic theory 
level, in the tracer limit one can assume that the velocity distribution function /(r, v;i) of the granular gas (excess 
component) obeys the (closed) Boltzmann equation while the velocity distribution function /o(r,v;t) of the tracer 
particles satisfies a (linear) Boltzmann-Lorentz equation. The problem is formally equivalent to consider an impurity 
or intruder immersed in a dilute granular gas, and this will be the terminology used in this paper. 

We consider a physical situation where the system (granular gas plus intruder) is in a steady state where weak spatial 
gradients of concentration, pressure and temperature coexist with a strong shear rate, which for the steady USF means 
strong dissipation [29(. Under these conditions, the resulting diffusion of intruder is anisotropic and, thus, tensorial 
quantities (Djj, D v , u ; and Dx^j) are required to describe mass transport instead of the conventional scalar transport 
coefficients [l7], HH, HH . Explicit expressions for the diffusion tensors Dij , D p ij and Dx,ij have been recently obtained 
[3(i| by solving the Boltzmann-Lorentz equation corresponding to the tracer particles by means of a perturbation 
expansion around a nonequilibrium sheared state (3ll - [33l | rather than the (local) equilibrium distribution j34|. This 
is the main new feature of this expansion (in contrast to the usual Chapman-Enskog method) since the reference 
state retains all the hydrodynamic orders (NS, Burnet, super-Burnett, • • • ) in the shear rate. As a consequence, the 
different approximations of this expansion are nonlinear functions of the coefficients of restitution as well as of the 
parameters of the mixture (masses and sizes) . 

The knowledge of the diffusion tensors allows us to study segregation by thermal diffusion. On the other hand, due 
to the anisotropy induced by the shear flow, a thermal diffusion tensor A is also required to characterize segregation 
in the different directions. Since in this paper we are interested in steady state conditions, wc consider a situation 
where the temperature gradient is orthogonal to the shear flow plane and parallel to gravity (i.e., d x T = d y T = 0, 
d z T 7^ and d x U y = a = const.). In this case, the segregation criterion is obtained from the thermal diffusion factor 
As, which is given in terms of the generalized mass transport coefficients D zz , D p ^ zz and Dt, zz - The use of these 
generalized non- Newtonian coefficients results in different and more general segregation criteria than those of previous 
works in sheared systems [2!|, which are limited to nearly elastic particles. 

Segregation is induced and sustained by both small gravity field and/or temperature gradient. The signature of 
A z provides a segregation criterion that shows a transition between upwards and downwards segregation (or BNE 
and RBNE when the intruder is large) by varying the parameters of the system. In particular, wc find that the form 
for the upwards/downwards segregation transition depends very sensitively on the value of gravity relative to the 
thermal gradient in such a way that depending on this one or the other mechanism vastly predominates in the space 
parameter. Moreover, our results show differences with those derived [I?], EH when the gas is driven by a stochastic 
thermostat that mimics the effect of a thermal bath J35|. These differences lead to interesting results, most notably, 
we found the upwards/downwards segregation mechanisms are much stronger when shear is input in the system. This 
may of course be a signature of more effective segregation process, with a direct impact for applications. 

The plan of the paper is as follows. First, the thermal diffusion factor A z is defined and evaluated in section ILTl by 
using a hydrodynamic description. In section ITO1 we determine the magnitudes needed to calculate A 2 (stress tensors 
of gas and intruder and zz elements of the diffusion tensors Dij, D p ^j and Dtaj)- All these quantities are explicitly 
obtained after solving the set of (inelastic) Boltzmann equations by means of the aforementioned Chapman-Enskog- 
like expansion. The knowledge of the above quantities yields A z as a function of the parameter space of the problem, 
namely, the mass and diameter ratios, the two independent coefficients of restitution for collisions among gas-gas and 
intruder-gas particles and the reduced gravity (gravity over thermal gradient). The form of the phase diagrams of 
segregation is investigated in section ITVl by varying the different parameters of the system. In addition, a comparison 
with the theoretical results [l3| derived when the system is thermalized by a stochastic thermostat is also carried out. 
Finally, we briefly discuss the results obtained in this paper in section |Vj 

II. HYDRODYNAMIC DESCRIPTION FOR THERMAL DIFFUSION UNDER SHEAR FLOW 

The model system considered is a dilute granular gas of smooth inelastic disks (d = 2) or spheres (d = 3) of mass 
m and diameter a, plus one intruder or impurity of mass tuq and diameter o~q. The presence of the intruder does 
not perturb the state of the granular gas and so, the model system (gas plus intruder) is formally equivalent to a 
dilute granular binary mixture in the tracer limit for the impurity component. Therefore, only gas-gas and intruder- 
gas particles collisions need to be taken into account. Collisions are inelastic and characterized by two independent 
(constant) coefficients of normal restitution a and a>o, respectively. We assume that the system (gas plus impurity) 
is in USF. This flow is characterized by constant densities n, n (number densities of gas and intruder respectively), 
uniform granular gas temperature T and a linear velocity profile U SjX = ay where a is the constant shear rate. In the 
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USF state the temperature changes in time due to the competition between two (opposite) mechanisms: on the one 
hand, viscous (shear) heating and, on the other hand, energy dissipation in collisions. A steady state is achieved when 
both mechanisms cancel each other and the fluid autonomously seeks the temperature at which the above balance 
occurs [29( | . 

The main goal of this paper is to study thermal diffusion of the intruder when the gas is under USF. We introduce 
small perturbations to our base state (USF). The perturbations are produced by weak gravitational field and small 
hydrodynamic gradients. These perturbations give rise to contributions to the mass flux, which can be characterized 
by generalized transport coefficients. Therefore, intruder segregation will be determined by the competition between 
these two different perturbations and by the relevant parameters values. 

The transport properties we need result from a general perturbation of the USF for which the flow velocity may be 
expressed as U = U s + SXJ where SXJ is a small perturbation. Here, XJ S = a • r with a = aSi X Sj y . Thus, under these 
conditions, the macroscopic balance equations for the system associated with this disturbed USF state are given by 

M 

d t n + XJ S • Vn + V • (nSXJ) = 0, (f ) 

d t n + XJ S • Vn = -V • (n SV) - (2) 

m 

dtSU + a • SXJ + (U s + SXJ) ■ VSXJ = -(mn)- 1 (V • P - nmg) , (3) 

i ndtT + ( L n {XJ s + SXJ) ■ VT + aP xy + V • q+ P : VSXJ = ~p(, (4) 

where we have assumed that the gravitational field g is the only external force. Here, g = —ge z , where g is a positive 
constant and e 2 is the unit vector in the positive direction of the z axis. Moreover, jo is the mass flux of intruder, P 
is the pressure tensor, q is the heat flux and ( is the cooling rate associated with the energy dissipation of collisions 
between gas particles themselves. It must noted that the balance equations (jT))— (J3|) can be exactly obtained from the 
(inelastic) Boltzmann equation and they provide the basis for developing a hydrodynamic description of this disturbed 
USF state. As usual, to get a hydrodynamic description one has to represent the fluxes jo, P and q as well as the 
cooling rate £ as explicit functionals of the hydrodynamic fields and their gradients. Once this constitutive equations 
are determined, the hydrodynamic equations ©-(H) become a closed set of equations for the fields n, n , XJ and T. 

Thermal diffusion is caused by the relative motion of the components of a mixture due to the presence of a thermal 
gradient. As a consequence of this motion, a steady state is achieved in which the separating effect arising from the 
thermal diffusion is balanced by the remixing effect of ordinary diffusion [36| . We intend to calculate segregation 
criteria for steady states that weakly deviate from the USF. Specifically, we are interested in a steady state with 
SXJ = and with (small) gradients of density, pressure and temperature only along the z axis. These small gradients 
are caused by boundary conditions for the temperature and a weak gravitational field. The hydrodynamic equations 
@-(U]) admit a steady solution for these states. 

Let us explain how we will obtain segregation criteria. From an experimental point of view, the amount of segregation 
parallel to the thermal gradient can be characterized by the thermal diffusion factor A z , which in a steady state is 
defined through the relation 



. 91nT d , 

A *^r- = ~7r ln 

az az 



(*)■ < 5 » 



We will consider also that the temperature gradient is directed downwards (dT/dz < 0); i.e. in the same direction as 
gravity. Thus, when A 2 > (i.e., d z \n(no/n) > 0), the intruder tends to rise with respect to the gas particles (BNE 
when the intruder is larger than the particles of the gas) and when A z < (i.e., d z \n(no/n) < 0), the intruder falls 
with respect to the gas particles (RBNE when the intruder is larger than the particles of the gas). 

We obtain now a relation of the type §5§ from the balance equations. In the steady state with SXJ = 0, the mass 
flux jo ;2 vanishes according to the balance equation ([2]). To first order in the spatial gradients, the mass flux j'o lZ is 
given by 



dxo m dp mn dT 

Jo, z = -moD zz — —D v . zz - — L> Tzz —, (6) 

az T az 2 az 



where Xq = n /n. Here, D ZZ1 D p zz and Dt.zz arG the zz elements of the diffusion tensor Dij, the pressure diffusion 
tensor D p jj and the thermal diffusion tensor Dx.ij, respectively (As we said, due to the presence of the shear flow, 
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diffusion process is anisotropic and thus tensorial quantities are required to describe mass transport under USF). 
Thus, the condition j Q z = yields 

dx _ m D PtZZ dp p D T , ZZ dT 

dz rnijT D zz dz m§T D zz dz 

Momentum balance equation @ for our steady state reduces to 

dP„, 



dz 



P9, (8) 



where p = ran is the mass density of the gas. This equation will allow to express the gradient of p as a function of 
the other two gradients and thus, to obtain a relation of the type ([5]) from the condition jo tZ = 0. 
In the hydrodynamic regime, the pressure tensor Pjj has the form [H, HH 

P l0 =pP* 3 (a*) 1 (9) 

where p = nT is the hydrostatic pressure and a* = a/v is the (reduced) shear rate. Here, 

_(d-l)/2 a 

-WJW^' (mT) (10) 

is an effective collision frequency. According to Eq. ([5]), the spatial dependence of P zz occurs explicitly through p and 
through its dependence on a*. As a consequence, 

dP** dp dP* zz da* 

dz zz dz 1 da* dz 

dv dT 1 v 

= £(\~ a *d a ,)P* zz + —-^a*{d a ,P: z X (11) 

where use has been made of the identity 

da* 71 



2 d z \nT-d z lnp) . (12) 



Using Eq. (jlip , we can obtain an expression of the gradient of p 

dz " P* z - a*{d a .P* z ) 
Use of (|llip into ([7]) and substitution of Eq. (J7} into Eq. ((SJ) finally leads to 



(13) 



A d t. zz - {PL - ^(da.p^y 1 d;, zz ( g * + \ a *{d a ,p* zz )) 

A - = ^ L , (14) 



where we have introduced the reduced coefficients 



D zz = -—Dzz, D pzz = 7^D p . zz , D Tzz = ——Dt.zz (15) 

and 

9* = JWr\ < ( 16 ) 

is a dimcnsionless parameter measuring the gravity relative to the thermal gradient. This quantity measures the 
competition between these two mechanisms {g and d z T) on segregation. 

The condition A z = provides the criterion for the upwards/downwards segregation transition. Since the diffusion 
coefficient D* zz is positive (as will be shown later), according to Eq. (fT4|) . the condition A z — implies 

(P; z - a*(d a *P: z )) D T ,zz = D* P ,zz {9* + \a{d a *P* zz )^ . (17) 

This equation delimits the upwards and downwards segregation regimes in a granular gas driven by shear flow. The 
explicit form of the pressure tensor and the diffusion coefficients relevant for this problem were calculated in previous 
works [1^, [3(1 IH, HH • We shortly explain in the next section the procedure followed to obtain them. 
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III. BOLTZMANN KINETIC THEORY 



We adopt now a kinetic theory point of view and start from the set of Boltzmann kinetic equations for the system 
(gas plus intruder). In this description, all the macroscopic properties of interest of the system are determined from 
the one-particle distribution function of the gas /(r, v;i) and the one-particle distribution function of the impurity 
/o(r, v; t). While the time evolution of the distribution function /(r, v; t) is given by the (closed) inelastic Boltzmann 
equation [37l. |38|. the distribution function /o(r, v;t) obeys the (linear) Boltzmann-Lorentz equation. 

In this paper we are interested in a flow that differs slightly from the steady USF [H, [2§|. Under these conditions, 
the kinetic equations for / and /o have been recently solved [30I433T ] by means of a Chapman-Enskog-like expansion 
around local shear flow distributions. Therefore, we look for solutions of the form 

/ = / (o) + / (D + ... ) / = / (o) +/ (D + ... ) (18) 

where the reference zeroth-order distribution functions /(°) and are the corresponding local versions of the 
USF distributions of the gas and the intruder, respectively. According to this perturbation scheme, the successive 
approximations /' fe ' and /g fe ' arc of order k in the spatial gradients of concentration xq, pressure p and temperature 
T but retain all the hydrodynamic orders in the shear rate a. As said in the Introduction, this is the main new 
ingredient of this expansion with respect to the conventional Chapman-Enskog method [34[ . This feature allows us to 
deploy for the first time a segregation theory based on non-Newtonian hydrodynamics for granular sheared systems. 

The rheological properties of the reference states /(°) and are related to the pressure tensors P^ ^ and P °\ 
respectively. They are defined as 

= J dvmVMfWiV), P ( °> = J dvmViVjf^iV), (19) 

where V = v — U s is the peculiar velocity. In addition, the steady state condition in the USF problem requires that 
the viscous heating term a*P* y is exactly compensated for by the collisional cooling term £* [28|, [2{|. In that case, 
the (reduced) shear rate a* and the coefficient of restitution a are not independent parameters but they arc coupled 
by the relation 

^Ky = ~C, (20) 

where P*j = pjp Jp and £* = Q/v. Equation (|20p shows the intrinsic connection between the shear field and dissipation 
in the system. In fact, in the elastic limit (a = 1), £* = and the (reduced) shear rate vanishes. The explicit 
expressions for the (reduced) pressure tensor P?-, Pq^ = P^/xop and their derivatives with respect to a* in the 
steady USF were obtained in Refs. [29|, HU, [33| . They are displayed in the Appendix [Al for the sake of completeness. 
Since a < 1 the range of (reduced) shear rates is defined, according to Eq. (|A3[) . in the interval < a* < (3 + 
2d)y/(d + 2)/(32d{d+l)). Thus, for hard disks (d = 2), < a* < 1.01 while < a* < 1.03 for hard spheres (d = 3). 
This clearly shows that non-Newtonian effects become important as the dissipation increases. To illustrate this non- 
Newtonian behavior, Fig. [T] shows a* versus a for hard disks and spheres. It is apparent that the (reduced) shear 



rate (which gives the steady granular temperature) is a nonlinear function of a in contrast to previous analysis [25 1 
carried out for nearly elastic particles where a* 2 oc 1 — a. 

To complement the results shown in Fig. [Tl the dependence of the relevant elements of P*j and Pq ^ on a is plotted 
in Fig. [2] for the system d = 3, <jq/<j = 2 and mo/m — 4. As expected, the transport properties in the steady USF 
state are inherently different from those of the NS description [which means here nearly elastic particles due to Eq. 
([20]) ]. A non-Newtonian signal of this behavior is the existence of the normal stress differences in the shear flow plane. 
It must remarked that, although the expressions of P*j and Pq ^ have been obtained by an approximate solution based 
on Grad's method, their a-dependence compare quite well with Monte Carlo simulations even for strong dissipation 

mil. 

Once the rheological properties of the gas and the intruder are well characterized, the diffusion coefficients D zz , 
Dp. Z z and Dt,zz can be obtained to first order in the expansion. The expressions of these generalized coefficients were 
derived in Ref. While the coefficient D zz decouples from the other two, the coefficients D p _ zz and Dt,zz obey a 
set of coupled equations. These equations are solved in the Appendix [B] The results show that D zz > while D P:ZZ 
and Dt, zz have not a definite signature. 

Before exploring the dependence of the parameter space on the form of the phase diagrams, it is instructive to 
consider certain limit situations. For example, there is no segregation when the intruder and the particles of gas are 
mechanically equivalent (mo = m, ctq = a, and ctQ = a). This is consistent with the results derived in the Appendix 
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FIG. 1: Plot of the reduced shear rate a* as a function of the coefficient of restitution a for d = 2 (dashed line) and d = 3 
(solid line). 




FIG. 2: Dependence of the diagonal elements (a) P^ x an d Po,xx, (b) Pyy and Po, vy , and the off-diagonal elements (c) P£ v and 
-Po.xi; on the (common) coefficient of restitution a = ao in the three-dimensional case for ao/a = 2 and mo/m = 4. The solid 
lines correspond to the elements PJ of the granular gas while the dashed lines are the elements Poij of the intruder particles. 



[B] since in this limit D PtZZ = Dt.zz = [the right hand side of eqs. (|B3|) and (|B4I) vanish since P* z = Pq zz ] and so, 
A z = for all values of a. Another reference situation is the elastic limit [a = ao = 1, which implies a* = in the 
steady state condition (|2T)|) ]. In this limit case, P zz = 1 and the diffusion coefficients behave as 



0, (21) 



where 



"° - flW)' 1 1,-2-- J V to (to + to ) (22) 
is a (positive) collision frequency. Consequently, for elastic collisions the segregation criterion (|17[) becomes 

9* 0. (23) 

V TO / 
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FIG. 3: Plot of thermal diffusion factor A z versus the (reduced) gravity \g*\ for for inelastic hard spheres (d — 3) in the case 
mo/m = ao/a = 2. Three different values of the common (a = ao) coefficient of restitution a have been considered: a = 0.9 
(solid line), a = 0.8 (dashed line), and a = 0.7 (dotted line). 
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FIG. 4: Plot of thermal diffusion factor A z versus common (a — cto) coefficient of restitution a for inelastic hard spheres 
(d = 3) in the case ao/a = f and mo/m = 2. Three different values of the reduced gravity \g*\ have been considered: \g*\ = f 
(solid line), \g*\ — 0.1 (dotted line), and \g*\ =0 (dashed line). 

In the absence of gravity, Eq. ([25)1 holds trivially, and so the intruder does not segregate. On the other hand, when 
\g*\ ^ the solution to (|23[) is simply mo = m, namely, segregation is predicted for particles that differ in mass, no 
matter what their diameters may be (l3| . For inelastic systems, in general the criterion (|17[) is rather complicated 
since it involves all the parameter space of the problem. 

Equation (|14p clearly shows that A z is a linear function of the reduced gravity \g*\. This is illustrated in figure 
[3] where thermal diffusion is plotted as a function of \g*\ for mo/m. = ao/a = 2 and several values of a common 
coefficient of restitution (a = ao). For the case represented here, downwards segregation is dominant, except for 
quite small values of \g*\ (i.e., there is a critical value \g*\ such that a transition upwards segregation =4> downwards 
segregation occurs for \g*\ > \g*\). It is to be noticed that the opposite behavior is present in the segregation results 
for dense systems driven by a stochastic thermostat [UHil- Figure 2] shows the cv-dependence of the thermal diffusion 
factor for different values of the reduced gravity \g*\. We observe that the presence of gravity changes dramatically 
the regions of positive and negative A z . In particular, for the case considered here and for not quite small values of 
a, we observe that the downwards segregation region increases with the reduced gravity. 

IV. PHASE DIAGRAMS FOR SEGREGATION 

As it can be noticed from Eq. (|14[) and from the Appendix A, the thermal diffusion factor A z depends on many 
parameters: the dimcnsionlcss gravity g* , the mass ratio mo/m, the size ratio o^ja, and the coefficients of restitution 
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FIG. 5: Phase diagram for segregation for inelastic hard spheres (d = 3) in the absence of gravity (\g*\ = 0). Three different 
values of the (common) coefficient of restitution a have been considered: a — 0.9 (solid line), a — 0.8 (dashed line), and a — 0.7 
(dotted line). Points above the curve correspond to A z > (upwards segregation) while points below the curve correspond to 
A z < (downwards segregation). 



a and arj. Given that the parameter space is fivefold, for simplicity, we take the simplest case of common coefficient 
of restitution a = an. This reduces the parameter space to four quantities. 

The zero contour of A z separates regions of positive (upwards segregation) and negative (downwards segregation) 
A z . Points lying on the zero contour correspond to values of the parameters of the system for which the intruder 
does not segregate. According to Eq. (|17[) . segregation is sustained by both gravity and the thermal gradient. The 
combined effect of both g and d z T on thermal diffusion is through the dimensionlcss gravity g* < defined by Eq. (|16l) . 
Although our previous perturbation analysis [3(| assumes that the external field is of the same order of magnitude as 
d z T, it is instructive to separate the influence of each one of the terms appearing in Eq. (|17p on segregation. Thus, 
two specific limit situations will be considered first in the next subsections: (a) absence of gravity (g = 0), and (b) 
homogeneous temperature (d z T = 0). 



A. Absence of gravity (|<?*| = 0) 

In this case, segregation of the intruder is due only to the presence of the thermal gradient. Under these conditions, 
| g*\ -> and Eq. ([P7]) reduces to 

P; z D* TiZZ = a*(d a ,P; z ) (d^ zz + \d* p>z )j . (24) 

Of course, this condition is trivially satisfied in the elastic case (for which a = 1, a* = 0, and D^, zz = 0). For inelastic 
systems (a ^ 1), the influence of each term in Eq. (|24]l is still intricate due to the presence of shear flow. As an 
illustration, figure [5] shows the phase diagram in the {mo/m, errj / 'cr}-planc for three different values of the coefficient of 
restitution a. As expected, when mo/m = (To/a = 1, the species are indistinguishable, and so A 2 = 0. Consequently, 
all zero contours of A z pass through the point (1, 1). It is apparent that, for en < c, the main effect of dissipation 
(or equivalently, the reduced shear rate a*) is to reduce the size of the upwards segregation region while the opposite 
happens when the intruder is larger than the particles of the gas. In the latter case, the influence of dissipation on 
the phase diagram is much smaller than in the former case (when <7o < a) so that all the curves tend to collapse in 
a common one for sufficiently large values of the size ratio cro/cr. We also observe that in the zero gravity limit large 
intruders will in general tend to move towards hotter regions, since the upwards segregation is dominant and occupies 
most of the parameter space: in fact, for large intruders (o'o/ cr > 1); upwards segregation is the only active mechanism 
for <tq/<t > 1.4. This is to be put in contrast with the case of no shear (vibrated systems in dense systems), where 
segregation tends to be predominantly of the downwards segregation type as the intruders get larger (l9j . just the 
opposite behavior seen here for sheared granular gases. 
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FIG. 6: (color online) Segregation phase diagram for inelastic hard spheres (d = 3) in the absence of thermal gradient 
(|<7*| — i> oo). Three different values of the (common) coefficient of restitution a have been considered: a = 0.9 (solid line), 
a = 0.8 (dashed line), and a = 0.7 (dotted line). Points above the curve correspond to A z < (downwards segregation) while 
points below the curve correspond to A z > (upwards segregation). The dashed-dotted line refers to the results obtained when 
the gas (for a = 0.8) is driven by a stochastic external force. 



B. Thermalized systems (d z T = 0) 

Let us consider now a system with both negligible temperature and mole fraction gradients (d z T 0, 8 z xq — > 0). 
Thus, the segregation of the intruders is only driven by gravity. This situation (gravity dominates over thermal 
gradient) can be achieved in vibrated or sheared systems in computer simulations and real experiments fjj, l43H45j . In 
this case (\g*\ — > oo), the sign of A z is the same as that of the pressure diffusion coefficient D* zz and so, the criterion 
(|17|) becomes simply 

D* p , zz = 0. (25) 

In the elastic case, D* zz oc m — too and so, the segregation criterion is too = to. For inelastic gases, the dependence 
of D* zz on the parameter space is quite complex. 

As we said in the Introduction, previous results for thermal diffusion have been obtained when the gas is driven 
by means of a stochastic external force (thermostat) that mimics the effect of a thermal bath. This external driving 
method is usually employed in computer simulations [5i| to compensate for cooling effects associated with the inelas- 
ticity of collisions. Under these conditions, an explicit expression for the thermal diffusion factor based on the NS 
transport coefficients has been recently obtained flo| . This expression for A z leads to the segregation criterion fl3j 

^ = ^, (26) 

1 TO 



where the temperature ratio is determined from the condition £o21)/ m o = QT/m. The segregation criterion (|26[) 
compares well with molecular dynamics simulation results for the case of the steady state of an open vibrated granular 
system in the absence of macroscopic fluxes [l2[ . As expected, it is clear that the segregation criterion ([25)1 (obtained 
from the NS description) differs from the one derived here for sheared gases (see Appendix |B|) . The discrepancies 
between ([25]) and ([26]) are a direct consequence of the inherent non-Newtonian features of the USF state not present 
at the NS level. 

A typical phase diagram for thermalized sheared systems delineating the regimes between upwards and downwards 
segregation is plotted in figure [SI The results obtained from the relation Ip26"]l (stochastic driving thermostat) for 
a = 0.8 are also shown for comparison. In spite of the differences between the criteria (|25[) and ([26]) . the shape of the 
phase diagram obtained from both relations agrees very well for large size ratios. In particular, when ao/a < 1, the 
main effect of inelasticity is to reduce the size of the downwards segregation region, while the opposite happens when 
<7q/(7 > 1. On the other hand, the effect of dissipation on the form of the phase diagrams is much more important when 
cto < (J than when uq > a. Moreover, comparison between figures [5] and [6] clearly shows that the presence of gravity 
changes dramatically the form of the phase diagram since the regions of positive and negative A z are interchanged. 

In the same way as it happens for \g*\ = 0, there is one clearly predominant segregation mechanism, the downwards 
segregation this time, being the only active mechanism for toq/to > 1 if a^ja > 1. Interestingly, the downwards 
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FIG. 7: (color online) Segregation phase diagram for inelastic hard spheres (d = 3) in the case \g*\ = 1. Three different values 
of the (common) coefficient of restitution a have been considered: a = 0.9 (solid line), a — 0.8 (dashed line), and a = 0.7 
(dotted line). Points above the curve correspond to A 2 < (downwards segregation) while points below the curve correspond 
to A z > (upwards segregation). The dashed-dotted line refers to the results obtained when the gas (for a = 0.8) is driven by 
a stochastic external force. 
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FIG. 8: Plot of thermal diffusion factor A 2 versus the coefficient of restitution ceo for inelastic hard spheres (d = 3) for the 
system ao/cr — mo/m = 1 when \g*\ — 1. Three different values of the coefficient of restitution a have been considered: a — 0.9 
(solid line), a = 0.8 (dashed line), and a = 0.7 (dotted line). 

segregation is not present for small but heavy intruders if the granular gas is sufficiently inelastic. This last situation 
is in contrast with the vibrated systems since there the downwards segregation still occurs in this region [l9j . 

C. General case 

We analyze the form of the phase diagrams for finite values of the reduced gravity \g*\. Figure [7] shows the phase 
diagram when \g*\ = 1 (gravity comparable to the thermal gradient) for the same cases as considered in the previous 
figures. In general, the form of the phase diagram for \g*\ = 1 is quite similar to the one obtained when \g*\ = oo; i.e., 
\g*\ = 1 is already sufficiently large to show the high gravity segregation behaviour (similar to that of the thermalized 
systems ) . Comparison with the stochastic driving results shows now better agreement with the sheared results than 
those observed for thermalized systems (see figure • 

All previous graphs have been obtained by assuming a common coefficient of restitution. Therefore, it is worth 
studying the (pure) effect on inelasticity on segregation in shear flows. In order to illustrate this case, we consider 
the case mo = m, and gq = a with \g*\ = 1, corresponding to a system of particles that differ by their coefficients 
of restitution only. Clearly, when all the coefficients of restitution are equal (a = ao), the thermal diffusion factor 
vanishes (A z = 0) and so, as expected the intruder does not segregate. Figure |S] presents plots of A z versus ao for 



\ 




Upwards segregation 




11 



different values of a. When the coefficients of restitution are different from each other, thermal diffusion is different 
from zero so that segregation in the presence of a temperature gradient can therefore occur as a consequence of 
inelasticity only. This effect is also present when the mixtures are vertically vibrated (Til l47l]. But, again contrary to 
what happens in the vibrated systems, in the sheared system the upwards segregation region increases for increasing 
inelasticity. 

Therefore, we have seen that in most cases the qualitative behavior of the segregation criteria for the sheared system 
is rather different from the criteria for vibrated systems. Moreover the shearing has the effect of clearly favoring just 
one segregation mechanism: upwards segregation for \g\ < \g c \ and downwards segregation for \g\ > \g c \. 



V. DISCUSSION 



The problem of segregation of an intruder in a sheared granular gas has been addressed in this paper. The relative 
motion of the intruder with respect to the particles of the gas is caused by the combined effect of gravity and a 
temperature gradient. In this case, thermal diffusion (or thermophoresis in its single-particle manifestation) forces 
large or massive particles to move down temperature gradients [481 ] . Under these conditions, the amount of segregation 
parallel to the thermal gradient can be measured by the thermal diffusion factor. However, the analysis of thermal 
diffusion in a strongly shearing gas is an intricate problem basically due to the anisotropy induced in the system by 
the shear flow. For this reason, in general a thermal diffusion tensor is needed to describe the segregation process in 
the different directions. In this paper, for the sake of simplicity we have assumed that the temperature gradient is 
orthogonal to the shear flow plane. 

Under the above conditions, the thermal diffusion factor A z [defined by Eq. ([5])] has been obtained in a steady state 
where weak spatial gradients of concentration, pressure and temperature directed along the vertical direction (parallel 
to gravity) coexist with a strong (constant) shear rate a. Our approach to determine the thermal diffusion factor 
follows two complementary routes. First, by using a hydrodynamic description, A z has been expressed in terms of 
the pressure tensor Pjj of the gas and the transport coefficients D zz , D pzz , and Dt, zz associated with the mass flux 
of the intruder. Second, we have adopted a kinetic theory point of view and have considered some previous results 
[30L l3lj obtained by solving the Boltzmann kinetic equation for the gas and the Boltzmann-Lorentz equation for the 
intruder by means of a Chapman- Enskog-like type of expansion. This allows us to compute the thermal diffusion 
factor as a function of the mass and size ratios, the coefficients of restitution for gas-gas and intruder-gas collisions 
and the reduced gravity g* = pg/nd z T < 0. Once the explicit form of A z is known, the condition A z = provides the 
segregation criterion for the transition upwards segregation (regions of positive A Z )<J4> downwards segregation (regions 
of negative A z ). This criterion is given by equation (|17|) in terms of the (reduced) zz-element of the pressure tensor 
P* z and the diffusion transport coefficients D* z , D* zz , and zz - In order to make the paper self-contained, the 
complete expressions for all the above quantities have been displayed in the Appendices E] and [B] 

Some previous works pjL, 25] on thermal diffusion segregation in sheared systems have been based on the elastic NS 
transport coefficients. These coefficients arc obtained by the Chapman-Enskog expansion around the local equilibrium 
distribution. Thus, their expressions are limited to nearly elastic systems in the USF problem. The present study takes 
the (local) shear flow distributions for the gas particles and the intruder as reference states in the Chapman-Enskog 
method so that, the corresponding transport coefficients derived from this expansion retain all the hydrodynamic 
orders (NS, Burnett, super-Burnett, • • • ) in the shear rate. As a consequence, since the coefficients of restitution 
and the shear rate are not independent parameters in the steady USF [see Eqs. ([2U]) and/or (|A3[) ]. our segregation 
criterion goes beyond the weak-dissipation limit. This is perhaps the main new added value with respect to previous 
theoretical results [ij], Hj|- Moreover, our theory also takes into account the influence of both thermal gradients and 
gravity through the reduced gravity g* . 

In order to illustrate the form of the phase diagrams in the mass and size ratio plane, a (common) coefficient 
of restitution a = ao has been assumed. Two specific situations have been mainly studied: \g*\ =0 (absence of 
gravity) and \g*\ — > 00 (homogeneous temperature). The shear field has the neat effect of selecting more strongly 
the upwards segregation mechanism: we have found that upwards segregation predominates very strongly for small 
gravity \g*\ (is in fact the only active mechanism for large and heavy intruders if gravity is small). These results are 
in contrast with segregation in the absence of shear, where there is no such a strong predominance of the upwards 
segregation mechanism for small gravities fl9| . Regarding the impact of the inelasticity of collisions on segregation, 
the results show that the influence of dissipation on thermal diffusion is more important when the thermal gradient 
dominates over gravity (g* = 0) than in the opposite limit (\g*\ — > 00), especially when the intruder is smaller than 
the particles of the gas. In addition, in general the effect of inelasticity on segregation mechanisms is less sensitive 
when the diameter of the intruder is greater than the particles of the gas. We also observe that the influence of gravity 
on the form of phase diagrams is quite important: while upwards segregation is dominant for small <?*, downwards 
segregation dominates if g* is large. To complement the previous study, we have also considered the case of identical 
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mass and size for the intruder and gas particles but different coefficients of restitution (a ^ o<o). As it happens for 
vibrated granular mixtures [Til [4tJ , the intruder can segregate in a sheared granular gas on the basis of differences in 
their inelastic properties. 

The results derived here for thermal diffusion have been obtained in the tracer limit xq — > 0. This limit precludes 
the possibility of analyzing the influence of composition on the thermal diffusion factor A 2 , where previous results 
for dilute gases (l3l | have shown that the effect of xq on A z can be significant in many situations. The study on the 
dependence of thermal diffusion on composition is an interesting open problem, ft is also apparent that the results 
presented here are relevant to make a comparison with numerical simulations and/or experiments. In this context, 
we hope that this paper stimulates the performance of such simulations/experiments to check the relevance of kinetic 
theory to describe thermal diffusion segregation under shear flow. We plan to work along the above lines in the near 
future. 
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Appendix A: Pressure tensor of the gas and of the impurity in the steady USF 



The non-zero elements of the (reduced) pressure tensor of the gas P*j = Pij/p arc given by [2 



p yy - p zz - ■ ■ ■ - p dd - 2d + 3 - 3 J ' P '* x ~ d ( d ty P vv> ( A1 ) 



P 



■nj 



Ad, d+ l +{d - 1)a a* 
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(A2) 



where the relationship between the reduced shear rate a* = ajv [where v is defined by Eq. (|ll)p] and the coefficient 
of restitution a is 



a* (a) 
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Moreover, the (reduced) cooling rate C* = C/ u is 



(A4) 



The derivative of P* z with respect to a* is [30|, |32J, [32 
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where 
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(A6) 



and A = {dP* y /da*^ is the real root of the cubic equation 

d 2 

2a* 4 A 3 + Ada* 2 (C + f3)A 2 + — (7C* 2 + U(*f3 + 4/3 2 )A + d 2 /3(C* + Py 2 (2(3 2 - 2(* 2 - PC) = 0. (A7) 

Apart from the pressure tensor of the gas, another relevant transport property is the partial pressure tensor Po 
of the impurity. This quantity along with its derivative with respect to a* is also needed to compute the diffusion 
coefficients. The elements of the reduced pressure tensor Pq^ = Po,ij/n a T are [3lj 
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where 7 = To /T is the temperature ratio and 
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Here, a = (a + (To)/2, = moT/mTo is the mean-square velocity of the gas particles relative to that of the impurity, 
M = mj (m + m ) and M = 1 — M = m / (m + Too). The temperature ratio 7 is determined from the condition 
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where the "cooling rate" Q = Co/v f° r the impurity is given by 
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The derivative of the elements of P * y with respect to a* is more involved. They can be written as [3C 
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Moreover, Y' = (dY/d-y), X' = {dX/d-y), X' a = {dY/d-f), (' Q = (dQ/dj) and dj/da* = A 1 /A 2 where 

Ai = d Qa* X ~ *o) { (\a* X - X J [xi - ^o%) - (x^o% - Po, yy + XA) } 

+2a* 2 [xPlyy + X(dP; y /da*)] , (A22) 
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Appendix B: Diffusion coefficients Dz Z , D* Z! , and zz 



In this Appendix we give the explicit expressions of the tensors D* z , D* zz and zz . The diffusion coefficient D* 
is given by [3(1 HH 
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In the region of parameter space explored, fi* 2 > so that the diffusion coefficient D* zz is always positive. 

The pressure diffusion D* zz and the thermal diffusion zz coefficients are coupled. They obey the set of algebraic 
equations [30| 
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where A is the real root of (| A7|) and the expressions of the nonzero elements of the pressure tensors P^* and P *jj 
and their derivatives with respect to a* are given in the Appendix |Al The solution to the set of equations (|B3|) and 
(|B4[) is elementary and provides the explicit expressions of zz and D* z in terms of m /m, cro/o", a and «o- The 



results show that while D* z > 0, the coefficients D* zz and zz have not a definite sign 
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